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Abstract. We consider the number of nodes in the levels of unlabelled rooted random trees 
and show that the stochastic process given by the properly scaled level sizes weakly converges 
to the local time of a standard Brownian excursion. Furthermore we compute the average and 
the distribution of the height of such trees. These results extend existing results for conditioned 
Galton- Watson trees and forests to the case of unlabelled rooted trees and show that they behave 
in this respect essentially like a conditioned Galton- Watson process. 



1. Introduction 

We consider the profile and height of unlabelled rooted random trees. This kind of trees is 
also called Polya trees, because the enumeration theory developed by Polya allows an analytical 
treatment of this class of trees by means of generating functions (see [43] ) . The profile of a rooted 
tree T is defined as follows. First we define the fc-th level of T to be the set of all nodes having 
distance k from the root (where we use the usual shortest path graph metric). Let Lk{T) denote 
the number of nodes of the fc-th level. The profile of T is the sequence {Lk{T))k>o- For a random 
tree this sequence becomes a stochastic process. 

The first investigations of the profile of random trees seem to go back to Stepanov [45] who 
derived explicit formulas for the distribution of the size of one level. Further papers deal mainly 
with simply generated trees as defined by Meir and Moon Note that simply generated trees 
are defined by the functional equation 

(1) y{x) = x(j){y{x)) 

for their generating function but can also be viewed as family trees of a Galton- Watson process 
conditioned on the total progeny. Kolchin (see [32l [33]) related the level size distributions to 
distributions occurring in particle allocation schemes. Later Takacs [46] derived another expression 
for the level sizes by means of generating functions. Aldous [1] conjectured two functional limit 
theorems for the profile in two different ranges which were proved in [131 123j . The first author 
[To] studied restrictions of the profile to nodes of fixed degree. An extension to random forests of 
simply generated trees is given by the second author [24] . 

Later other tree classes have been considered as well. The profile of random binary search trees 
has been first studied by Chauvin et al. [5^ and later by Drmota and Hwang (see JJj and 
Random recursive trees have been investigated recently by Drmota and Hwang [17] and van der 
Hofstad et al. [48]. Related research was done by Chauvin et al. [5], Fuchs et al. [22], Hwang 
[29] [30] , Louchard et al. [34] , and Nicodeme [39] . Extremal studies of the profile (called the width 
of trees) of simply generated trees have been started by Odlyzko and Wilf [40]. The distribution 
including moment convergence has been presented independently in Chassaing et al. [J and the 
authors of this paper [15] . For other tree classes we refer to the work of Devroye and Hwang [9] 
and Drmota and Hwang [17) . A general overview on random trees which also strongly highlights 
the profile of trees can be found in the first author's book on random trees [12j . 
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Whereas simply generated trees have an average height of order y^, the other tree classes 
mentioned above have height of order logn. Polya trees do not belong to the class of simply 
generated trees. Since we are not aware of any rigorous proof of this assertion in the literature, 
we will present a (rather simple) proof of this in the next section. To our knowledge, so far the 
fact was only underpinned by the following argument concerning the generating functions of both 
tree classes. The argument works for many tree classes considered in the literature, e.g. Cayley 
trees, plane trees, Motzkin trees, binary trees, and many more. In all these cases the function (/)(?/) 
in ([T]) is entire function or meromorphic. But it is not at all clear that this argument is still true 
if (f>{y) has a more complicated singularity structure. In these cases (i.e., entire or meromorphic 
4>{y))i the generating functions enumerating the number of simply generated trees and Polya trees, 
respectively, have a fundamentally different singularity structure. Whereas the first one has one or 
a finite number of singularities (the latter occurs in the periodic case) on the circle of convergence 
and allows analytic continuation to a slit plane (with the possible exception of finitely many 
isolated singularities which are of algebraic type even if the function itself is not algebraic), the 
generating function associated to Polya trees is much more complicated. In fact, for the latter 
function the unit circle is a natural boundary (i.e., no analytic continuation beyond it is possible). 
There is exactly one singularity on the circle of convergence of the power series expansion at 0, but 
the analytic continuation has an infinite number of singularities inside the unit circle. Each point 
on the unit circle is an accumulation point of the set of singularities. These facts follow from the 
functional equation defining this generating function and the fact that the power series expansion 
around zero has radius of convergence strictly smaller than one (see next section) . It also involves 
an analytically complicated structure like the cycle index of the symmetric group. Due to this 
difference with respect to the analytic behaviour of the generating function Polya trees are not 
simply generated and therefore they cannot be represented as branching processes. 

Note that the rather complicated singularity structure does not affect the asymptotics of sta- 
tistical parameters like the number of trees, the profile or the height. The behaviour of these 
parameters is determined by the dominant singularity which is, as all other singularities as well, 
an isolated singularity, exactly as in the case of simply generated trees. Indeed, Polya trees behave 
in many respects similar to simply generated trees (compare with [331 [571 133 ISHl [131 [IS]) or the 
recent work of Marckert and Miermont [35 who showed that binary unlabelled trees converge in 
some sense to the continuum random tree, i.e., the same limit as that of simply generated trees. 
Moreover, Broutin and Flajolet [3 showed -^/n-behaviour for the height of binary unlabelled trees. 
Hence it is expected that the order of the height is ^/n as well. In this paper we will give an 
aflirmative answer to this question. This justifies the choice of ^Jn for the scaling of the level sizes 
in the subsequent theorems. 

The plan of the paper is as follows. In the next section we present our main results. Then we 
will set up the generating functions for our counting problem of trees with nodes in certain levels 
marked. This function is given as solution of a recurrence relation which has to be analyzed in 
detail. Knowing the singular behaviour of the considered generating functions allows us to show 
that the finite dimensional distributions (fdd's) of the profile, i.e., the distributions of the sizes of 
several levels considered simultaneously, converge to the fdd's of Brownian excursion local time. 
The singularity analysis is carried out in Section [4] and the computation of the fdd's in Section [5l 
In order to complete the functional limit theorem we need to prove tightness. This means, roughly 
speaking, that the sample paths of the process do not have too strong fiuctuations (see ^ for the 
general theory). 

In the final section we turn to the height. The pioneering work on this topic was done by 
Flajolet and Odlyzko "20^ and Flajolet et al. [19] in their studies of simply generated trees where 
they completed the program started in [8]. What we have to do is to show that the generating 
function appearing in the analysis of the height has a local structure which is amenable to the steps 
carried out in [30] and [TH]. This is done in the last section and leads to average and distribution 
of the height. 
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2. Preliminaries and Results 

First wc collect some results for unlabelled unrooted trees. Let denote the set of unlabelled 
rooted trees consisting of n vertices and y„ be the cardinality of this set. Polya [43j already 
discussed the generating function 

TI>1 

and showed that the radius of convergence p satisfies < p < 1 and that x — p is the only 
singularity on the circle of convergence l^j = p. He also showed that y{x) satisfies the functional 
equation 



(2) y{x) — X exp 



y{x') 



i>i 



Nowadays, this functional equation is easily derived by using the theory of combinatorial construc- 
tions which is presented in the comprehensive book of Flajolet and Sedgewick [21 . Indeed, a Polya 
tree can be viewed as a root with a multiset of Polya trees attached to it. Then the functional 
equation ([2]) pops out immediately from the multiset construction and its generating function. 
This functional equation can be used to compute the coefficients: 

(3) y{x) =x + x^ + 2x^ + \x^ + ^x^ + 20a;^ + 48a;^ + 115a;^ + 2862:^ + 719a;^° + . . . 
Later Otter flT showed that y{p) — I as well as the asymptotic expansion 

(4) yix) ^l-b{p- xfl'^ + c{p - + d{p - xf^ + ■■■ 
which he used to deduce that 

(5) y„ - ^^n-'l^p-. 

Furthermore he calculated the first constants appearing in this expansion: p 0.3383219, b ~ 
2.6811266, and c = /i w 2.3961466. 

We will return to the function y{x) in Section|4]and list a couple of useful properties in Lemma[l] 
after introducing some notations. 

Theorem 1. Polya trees are not simply generated. 

Proof: Let us assume that Polya tree are simply generated. Then the generating function y(x) given 
by must have a representation in the form ([l} where is a power series with non-negative 
coefficients. By ([3]) the functional inverse y~^{x) exists and we have y{x) ~ x and y~^{x) ~ x, as 
X — 0. This implies (/)(0) — 1. Plugging y~^(a;) into ([T]) we obtain 

X = y{y~^{x)) = y^^{x)(j){y{y~'^{x))) = y^^{x)(j){x) 

and consequently 

X X 



i.x) = — 



X ~ x'^ + x'^ - x^ + x^ - Ax"^ + llx^ - ISx^ -I- 18x10 -|- . . . 
^l + x + x^ - x^ -Vix^ -bx'' + 7x^ ~8x^ + x^° + ... 

which violates the requirement of non- negative coefficients for (t>{x). □ 



Remark. The sequence of the coefficients oiy{x) is A000081, that of y ^(x) is A050395 in Sloane's 
On-line Encyclopedia of Integer Sequences [?] 

The height of a tree is the maximal number of edges on a path from the root to another vertex 
of the tree. It turns out that the average height is of order i/n. 
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Theorem 2. Let denote the height of an unlahelled rooted random tree with n vertices. Then 
we have 

(6) Ei/„^T^x/^. 



and 

(7) Ei?; ^ (^^y r(r- l)r(r/2)C(r)n'-/2 
for every integer r > 2. 

The proof of this theorem is deferred to the last section, since the proofs of the auxiliary lemmas 
which will eventually establish the assertion will utilize similar techniques as needed to prove the 
next three theorems. 

Remark. Note that more information on the limiting distribution is available. Indeed, a local limit 
theorem holds as well. Let ylj^"^ denote the number of unlabelled rooted trees with n vertices and 
height equal to h and let 5 > arbitrary but fixed. If we set (3 — ^^fnjhh^, then, as n — )• oo, we 
have 

IJn \ n 

m>l 

uniformly for ^^j^g ^ < < S\/log n. A rigorous proof of this theorem was given by Broutin and 
Flajolet [3] for binary unlabelled trees. They also provide a moment convergence theorem, a weak 
limit theorem as well as large deviation results. 

Let Ln{t) denote the number of nodes at distance t from the root of a randomly chosen unla- 
belled rooted tree of size n. If t is not an integer, then define Ln{t) by linear interpolation: 

(8) L„(i) = (LiJ+l-<)L„(W) + (t- t>0. 
We will show the following theorem. 

Theorem 3. Let 

ln{t) = —J=Ln (t^/n) . 

Then ln{t) satisfies the following functional limit theorem: 

^(^^'(^■' _ 

in C[0,oo), as n oo. Here b and p are the constants of Equation Q and l{t) denote the local 
time of a standard scaled Brownian excursion. 

In order to prove this result we have to show the following two theorems 

Theorem 4. Let b, p, and ln{t) be as in Theorem\^ then for any d and any choice of fixed 
numbers ti, . . . ,td the following limit theorem holds: 

as n — > oo. 

Theorem 5. For all non-negative integers n, r, h we have 

(9) ^[Ln{r)-Ln{r + h)f <Ch'n 
where C denotes some fixed positive constant. Consequently, the process ln{t) is tight. 
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3. Combinatorial Setup 

In order to compute the distribution of the number of nodes in some given levels in a tree of 
size Ti we have to calculate the number yk\ra\k2m2---kci7ndn 

of trees of size n with rrii nodes in level 

fci , i = 1 , . . . , d and normalize by j/„ . 

Therefore we introduce the generating functions yk{x,u) defined by the recurrence relation 

yo{x,u) = uy{x) 

(10) yfc+i(a;,ti)=xexp|^|]^^^^^j , fc > 0. 

The function yk{x,u) represents trees where the nodes in level k are marked (and counted by 
u). If we want to look at two levels at once, say k and ^, then we have to take trees with height at 
most k and substitute the leaves in level k by trees with all nodes at level t — k marked (counted 
by v) and marking their roots as well (counted by u). This leads to the generating function 

yk/{x,u,v) = yk,e-k{x,u,v) satisfying the recurrence relation 

yo^i{x,u,v) = uyi{x,v) 

(11) yk+i,e{x,u,v) = xexp \ ■. , fc > 0. 

V.>i ' J 

In general we get therefore 

yki,...,ka{x,Ui, . . . ,U(i) ^ ^ ykimik2m2---kaman'iJ-i ' ' ' '''^ 

■mi,.. . ,md,n>0 

= yki,k2-kiM-k2,...,ki-kd-lix^ Ui,.. .,Ud) 

where 

yo,e2,-,idi^^ ui, • • ■ , Ud) = uiyi^^,„^i^{x, U2, • • • , "d) 

?/fc+i,^2,...A(a^:"i: ■■■,Ud) = a^exp I -. ^ j , fc > 0. 

The coefficients of these function are related to the process L„(t) (see ([5]) and the lines before) by 
P {i„(fci) = mi, i„(A:i + k2) = m2, . . . , i„(A;i + k^ + ■ ■ ■ + kd) ^ ma} = 

Vn 

where the ki are integers and the probability space of the measure P is the set of Polya trees with 
n vertices equipped with the uniform distribution. 

As claimed in Theorem[3j the process ln{t) = -^Ln{t) converges weakly to Brownian excursion 
local time. From f28^ (cf. f6', ITS] as well) we know that the characteristic function 4>{t) of the total 
local time of a standard Brownian excursion at level k is 



(12) 0(t) = 1 + 



^/t^ J y/—s cxp{k^/ —2s ) — sinh (k-\/~2s ^ 



7 



where 7 = (c — ioo, c + ioo) with some arbitrary c < 0. The characteristic function of the joint 
distribution of the local time at several levels ki, . . . , Kd was computed in [13] (for d — 2 already 
in [6] albeit written down in a form which does not exhibit the recursive structure) and is given 

by 

(13) (ti,..., t<j) = 1 + -^ / /Ki....,«:^(a;,ii,...,td)e"^da;, 

where 

/ki ,...,Kp {Xjil, ■ ■ ■ ^td) ^ Kl{x,tl + '3/^2 — KI (■ ■ • ^Kd_i — Kjj_2 (-^J ^d— 1 + '^'k^ — Kd_i {x , t d)) ' ' ' )) 
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with 

it^—x exp(— K-y/— 2x ) 



-xexp(/tV— 2a;) — ity/2smh 2x) 

For further studies of this distribution see [26l |42j |47] . 

In order to show the weak hmit theorem we have to show pointwise convergence of the charac- 
teristic function (f)kx,--- ,ka,n{ti, . ■ . ,td) of the joint distribution of ^L„(A:i), . . . , -^Ln{kd) to the 
corresponding characteristic function of the local time in some interval containing zero. We have 

Vn ^ 

This coefficient will be calculated asymptotically by singularity analysis (see [TS]) for kj = 
\_Kj^/n\. Thus knowing the local behaviour of ykix,u) near its dominant singularity is the crucial 
step in proving Theorem [3l This is provided by the following theorem will be the crucial step of 
the proof. 

Theorem 6. Set Wk{x,u) = yk{x,u) — y{x). Let x = p(l + u = e**/^, and k = [Ky/n\. 
Moreover, assume that \ args| > 9 > and, as n — oo, we have s — O (log^ n) whereas t and k 
are fixed. Then Wk{x,u) admits the local representation 

iP'p it^/^ exp {— nby/ ~ ps ) 



(14) Wk{x,u) 



2V^ y3i_M(i_exp(-K&v^)) 



um 



iformly for k = O {\/n). 



The proof is deferred to the next section. 

Note that Theorem [6] implies Theorem |4] for the case d= 1. If we set d = 1 in Theorem IH then 
we have 



<^fc^„(t) = (x,e^'^^) 



yn 

1 /" / it/./7T\ dx 



yk (x,e'"^) 



1 + TT- — / Wk {x,e"' '" ' '^^ 



(x, e^*/^) 



27rij/„ "'^ V"' " ) a;"+i 
where F is a suitable closed contour encircling the origin. Using Theorem [5] it is easy to show that 
'I^K^-nif) converges to (l)^b^/(2^){tb^/2^/2) as desired. 

The higher dimensional case is more involved, but relies on the same principles. A complete 
proof is given in Section [5l 

4. The Local Behaviour of yk{x,u) - Proof of Theorem [H] 

4.1. Notation. We will provide some frequently used notations now. 
We will study the local behaviour of yt by analyzing the quantity 

Wk{x,u) = yk{x,u) - y{x) 

which frequently involves the term 

Wk{x\u'') 



^k(x-u) := ^ ■ 



Furthermore, estimates of the partial derivatives 



7fc(a;,u) = —yk{x,u) and ll^\x,u) = —yk{x,u), i > 2, 



will be needed. 
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Figure 1. The domains A, A^, and Q. For analyticity arguments we will need 
the domain A. Sometimes, the proof goes without change for the larger domain 
as well. In that case we will use G for the sake of more generality. For arguments 
where the asymptotic behaviour near the singularity is important the domain A^ 
is used. The domains (not depicted here) are needed whenever the asymptotic 
behaviour for u ^ 1 is considered. Here uniformity is often important such that a 
universal e-neighbourhood of it = 1 which is independent of k does not serve our 
needs. 

The asymptotic analysis of Wk (resp. yk) enables us to apply Cauchy's integral formula and get 
the coefficients of yk{x,u) asymptotically (see the proof of Theorem 2] in the next section) which 
eventually leads to an integral of the form (|T2)) . Therefore estimates for y{x), provided in the next 
lemma, and the other functions appearing in our analysis are needed. The estimates will be valid 
in various domains. Therefore let us introduce 

(16) {x eC:\x\ < p + I arg(a; - p)\ > 0}, 

(17) A^ = {xeC:\x-p\<e, \ aig{x - p)\ > 6} 

Q = {x eC:\x\ < p + r], \ a.Tg{x - p)\ =^ 0} 

(18) Ek = {ueC:\u\<l,k\u-l\<fj} 
with e, 77, ?7 > and < 6* < |. 

Remark. In all the arguments which we will use in the following proofs it is always assumed that 
e and rj are sufficiently small even if it is not explicitly mentioned. 

4.2. Analysis of the local behaviour of yk{x,u). Obviously, Wk{x,l) = 0. Since yk{x,u) rep- 
resents the set of trees where the vertices of level k are marked, we expect that liuik^oo yk{x, u) — 
y{x) inside the domain of convergence. This is not obvious, but follows from what we derive in 
the sequel. We start with a useful property of y{x). 

Lemma 1. Provided that rj in (|16p is sufficiently small, the generating function y{x) has the 
following properties: 

a) For x € A we have |y(a;)| < 1. Equality holds only for x = p. 

b) Let X ^ p (1 — o,nd \t\ < Clog n for some fixed C > 0. Then there is a c > such 
that 

I ^ M ^ 1 / max(l,|t|) 
\y[x)\ < 1 - cy . 

c) For \x\ < p we have \y{x)\ < y{\x\) < 1. Moreover, near x = the asymptotic relation 
y{x) ~ X holds. 

d) There exists an e > Q such that 

(19) |y(x)|>min(^|,M^ 
for all X ^ 8. 
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Proof: The first statement, when restricted to \x\ < p, follows from the facts that y{x) has only 
positive coefficients (except yo = 0), y{p) = 1 and there are no periodicities. Extension to A is 
easily established by using (|4]) and continuity arguments. 

The second statement is an immediate consequence of the singular expansion (|4]) of y{x). 

The first inequality of the third statement follows from the positivity of the coefficients j/„. The 
same fact also implies that y{x) is strictly increasing in the interval [0, p] and therefore bounded 
by y{p) = 1. The asymptotic relation near follows from yo = and yi = 1 or from the functional 
equation ^ 

Finally, for proving the last statement we split the circle |x| < p into the smaller circle < e 
(let us call it C) and the annulus A := {x : e < |x| < p}. The function y{x) is analytic in 
CUA except at a; = p, but still continuous there. Since the annulus A is compact, |2/(a;)| attains a 
minimum there. By the functional equation ^ this minimum must be positive, since both factor 
on the right-hand side of ^ are nonzero. Since y{x) ~ x, as x — 0, a sufficiently small e guarantees 
that min|2,|^g 12/(2^)1 > £/2. Since this inequality holds for every smaller e as well, we can choose 
an e such that 

£ \x\ 

min |y(a;)| > - and \yix)\ > — for |a:| < e 

e<|2;|<p 2 2 

hold. By continuity this can be extended to (|19p . □ 

Next we derive an a priori estimate of in a small domain. 

Lemma 2. Let \x\ < + e for sufficiently small £ and \u\ < 1. Then there exist a constant L 
with < L < 1 and a positive constant C such that 

\wkix,u)\<C\u-l\-\x\-L'' 

for all non-negative integers k. 



Proof: We first note that by using the recurrence relation (|T0|) we obtain 

Wk+i {x, u) = yk+1 {x, u) - y{x) 



;exp ^yk{x\u'')^ 



y{x) 



(20) =?/(x) exp u;fe(x,M) + ^ 



Wk{x\u'') 



i>2 



For fc = we have \wo{x,u)\ = |u — 1| • \y{x)\ < C\u — l\\x\ since y{x) — O {x) as x ^ 0. We will 
then use the trivial inequality 

(21) |e--l|< l"^' 



1 _ M 

^ 2 



for the induction steps. However, in order to apply this tool we need some a-priori estimates. 
Obviously we have for |a;| < p and |w| < 1 

\wkix,u)\ < 2y{\x\) 

and consequently 



E 

i>i 



Wk{x\u') 



< 



2V«^21og«. 



i>l 



Since the function is convex for < a; < p and y{p) — 1 we get ^^j^ < 1 + because of the 
value of p. Consequently 



x\ \ p / p 
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Thus, if |a;| < + e (for a sufficiently small e > we have 



By using ([21]) we thus obtain 



exp I ^ 

i>i 



E 



Wk{x\u') 



<2p + 2- 



P 





1 




< 


^ p 



E 

i>i 



Therefore, if we assume that we already know \wk{x, u)\ < C\u — l||a;|L'^ (for \x\ < p'^ + e, \u\ < 1 
and some L with < L < 1) then we also get 



\wk+iix,u)\ < \y{x)\ 

\y{^)\ 



exp 



E 



Wk{x\u') 



- 1 



< 



E 

i>l 



i>l 
Wkix\u') 



p i>l 



1-p- 



l - X 



By convexity we have y{x) < xj p for < a; < p and, thus, there exists e > with y{p^ + e) < p. 
Consequently we get for |a;| < + e the estimate 

\wk+x{x,u)\ < CL'L''\x\\u-1\ 



holds where 



L' 



l-p-^\{l-p'-e) 



The value of L' is smaller than 1 if e > is sufficiently small. Thus, an induction proof works for 
L = L'. ' □ 

Corollary 1. For |u| < 1 and \x\ < p + s (e > small enough) there is a positive constant C 
such that (for all k > 0) 

\^kix,u)\ < C\u-1\L'' 
with the constant L from the previous lemma. 

Proof: We have 



1- X 



< C\u-1\L'' 



C\u - 1|L* 



□ 



Corollary 2. Let u £ ^fc and a; e Ag. Then 

Y,^k{x\u^)^0{L^) 



i>2 
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Proof: The assertion follows immediately from the previous corollary and Taylor's theorem. □ 

We will eventually need precise upper and lower bounds of Wk{x,u) near its singularity. By 
Taylor's theorem Wk{x,u) can be expressed in terms of ^k{x,u). Hence, let us consider 7^(2;, u) 
now. 

Lemma 3. For x £ Q (where rj > is sufficiently small) the functions Jk{x) can be represented 
as 

(22) jk{x):=-/k{x,l) = Ck{x)y{x)\ 

where the Ck{x) form a sequence of analytic functions which converges uniformly to an analytic 
limit function C{x) (for x Cz O) with convergence rate 

(23) Ck{x)=C{x)+0{L''), 

for some L with < L < 1. Furthermore we have C{p) — \b^P- 
There exist constants e,6,fi > and < ^ such that 

(24) \jk{x,u)\=0{\y{xt+') 
uniformly for x G and m G Sj. . 

Remark. Note that the constants e, ^, 77 > determine the domains and S^. Thus the theorem 
guarantees that there are suitable domains A^ and S^. of the shape PT)) and ^TE\\ where the 
estimates are valid. 

Proof: The first statement we have to show is that the functions 7^(2;) are analytic functions in Q. 
We prove this by induction. Obviously, the assertion holds for k = since 70(2;) = y{x). Assume 
it is true for k. Then Tk{x) := J2i>2^k{x'^) is analytic for \x\ < y/p and hence also in Q. Using 
the recurrence relation of yk(x,u), Equation (|10l) . we get 

7fe+i(x,w) = Axe^'=(-^")+s.(-,-) 
on 

(25) = yA;+i(a;, u)^7fc(a;*,u')w*"\ 

i>i 

This implies 

lk+i{x) = y{x)jk{x) + y{x)rk{x) 
which finally implies that 7^+1 (x) is analytic in Q as well. 

By solving this recurrence we obtain also the analyticity of Ck{x) — Furthermore, we 

will show that the sequence {Ck{x))k>o has a uniform limit C'{x) which has the desired properties. 

Setting w = 1 we can rewrite (|25p to 

(26) Ck+i{x)y{x)'^+' = Ck{x)y{x)'^+' + y{x) {Ck{x^)y{x^f+^ + Ck{x^)y{x^f+^ + ...). 
resp. to 

(27) Ck+i{x) = y^Ck{xy'^ ^ 

Set 



V- |y(a;')l'=+' 



If 77 > is sufficiently small then due to Lemma [T] we have 

sup^^<l foralH>2and sup = 0(1*) 

xdB \v{x)\ ^ee \y[x)\ 
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for some L with < L < 1. Consequently we also get 



for some L with < L < 1 (actually we can choose L = L ). Thus, if we use the notation 

11/11 sup \J{x)\ then ^ yields 
see 

(28) 

and also 



\\Ck+i\\<\\Cd{l + Lk) 



(29) \\Ck+i-Ck\\<\\Ck\\Lk. 

But pS]) implies that the functions Ck{x) are uniformly bounded in the given domain by 

\\Ck\\ <co :=n(l + ^^)- 



e>i 



Furthermore, (j29p guarantees the existence of a limit limk^ao Ck{x) = C{x) which is analytic in 
Q; and we have uniform exponential convergence rate 



\Ck-C\\ <co5]l, = 0(L^) 



e>k 



Hence, we get as desired. 
Finally, note that (for < p) 



^ 7k{x, 1) = X! ^^2^"^" = xy'{x). 



k>0 



n>l 



On the other hand, 



fe>0 fe>o 
Since C{x) is continuous in Q this implies 



1 - y{x) 



J2 1) = Y.(^ix) + O {L^))y{xf+'n = ^^^^^ + O (1) . 



c{p) = hm -y'i-)(]-y(-)) ^ ^ 

x^p y[x) 2 

where we used (jH). 

Let us turn to the second assertion. In order to obtain the upper bound we set for t <k 



(30) 



Ct= sup \-^i{x,u)y{x) \. 



Observe that by Lemma [T] 1 7^(0;, u)| < Cgsnp^^^^ \y{x)\'' < Ce- Therefore, by Taylor's theorem 
and Corollary [1] we obtain 



\ye+ 



i{x,u)\ < |?/(a;)|exp ^ 



\wiix\u')\ 



(31) 

By ([25l) we have 



< |y(:E)|eC?^l"-i|+0(^'). 



Ce+i = sup 



yi+i{x,u) 



Applying ((3T|) and the estimate in Corollary [2] we get 

gCf|«-l|+0(L'=) 



Ce+1 < sup 
xeAc, tie 



(32) 



y{xy 

<e'^«/'=C,(l + 0(LO) 



-ft;{x,u) +^-fe{x\u')u' 1 
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Set 

co = l[{l + 0{U)) 

and choose 77 > such that e^'^"'' < 2. We also choose e < fj and Q < 8 < ^ such that \y{x)\ < 1 
for I2; — p| < £ and | arg(x — p)\ > 9. li k > is fixed, then it follows by induction that, provided 
that |u- 1| < v/k, 

Ce < Y[{1 + OiU)) ■ e^'^ofii/k < {I < k). 

This completes the proof of the lemma, since Cq = 1 . □ 

The representation (I22p from Lemma [3] gives us a first indication of the behaviour of Wk{x, u) 
for u close to 1. We expect that 

(33) Wk{x, u)k{u^ l)-ik{x, 1) - (m - l)C{x)v{xf. 

This actually holds (up to constants) in a proper range for u and x, although it is only partially 
true in the range of interest (see Theorem [6]). 

In order to make this more precise we derive estimates for the second derivatives 7^ {x,u). 

Lemma 4. Suppose that \x\ < p — rj for some 77 > and \u\ < 1. Then uniformly 

(34) ^^^\x,u) = 0{y{\x\f+'). 
There also exist constants e, 0, 77 > such that 

(35) j^\x,u) = 0{k\y{xt+^) 
uniformly for u G Sfe and x G . 

Proof: By definition we have the recurrence (compare with (1251) ) 

l^Ux,u) = y,+,{x,u)J2^J^^\x\u^)u^^-' 
i>i 



(36) +yk+i{x,u)iY,n{x\u'X-' 

\i>l 

+ yk+i{x, u) y^{i - l)7fe(a;*, w*)u 



i>2 



rn] 

with initial condition 7q (x) — 0. 

[2] [21 

First suppose that \x\ < p — rj for some 77 > and \u\ < 1. Then we have |7[. (x, u)\ < Tk {\x\, !)• 
Thus, in this case it is sufficient to consider non-negative real x < p — rj. We proceed by induction. 
Suppose that we already know that 7|f' (x) :— j^^\x,l) < Dky{x)'^^^ (where Dq = 0). Then we 
get from ([36]) and the already known bound 7fe(x, 1) < Cy{x)^ from Lemma [3] the upper bound 



itUix) < Dky{xY+' + Dky{xY+'y^i 



^ y(x)^+i 



+ y{x)\c^y{xf'^+Y.^k{x\l) 



+ Cy{xf+^Y.^^-l) 



y{x') 



k+l 



y{x) 



k+l 



i>2 

< y(xf+^ {Dk{l + OiL")) + C'yip - 77)^= + 0(i'=)) 
where we used Corollary [5] in the last step. Consequently we can set 

Dk+i = Dk{l + 0{L'')) + Cyip - Tif + 0{L^) 
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and obtain that = 0(1) as fc — > oo which proves ([M)) . 

Next choose e, 6, fj as in the proof of Lemma [3] and set (for £ < k) 



Di — sup 



Ti ix,u)y{xy 



By the same reasoning as in the proof of Lemma |31 where we use the already proved bound 
\^e{x,u) < C|y(a;)|^+^, we obtain 

Di+i < De e^^l^(\ + 0(L^)) + e^^'^ + O(L^), 

that is, we have 

with ai = e''<^/'=(l + 0(LO) and = e^^l^ + ©(L^. Hence we get 

fc-i fc-i 

Du<Y.^, n 

J=0 i=j + l 



< k max/3j e'''^ [|(1 + 0{L% 



£>0 

= 0(fc). 

This completes the proof of psp . □ 

Using the estimates for 7^(2;, m) and 7^ (x, u) we derive the following representations for Wk{x, u) 
and Sfc(a;, w). 

Lemma 5. Let 8,6,7] and Ck{x) — ^k{x, l)/y{x)^^^ as in Lemma\^ Then we have 

(37) wu{x,u) = Cu{x){u-l)y{xf+^{l + 0{k\u-l\)), 

uniformly for u d 'E.k and x € A^. 

Furthermore we have for \x\ < p + rj (for some rj > 0) and \u\ < 1 

(38) ^k{x, u) = Ck{x){u - l)2/(x2)^+i +0{\u- Ipydx^'^) , 
where the analytic functions Ck (x) are given by 



(39) Cfc(a:)=^Cfc( 



i>2 



^y{x^) 

They have a uniform limit C{x) with convergence rate 

Ck{x)^C{x)+0{L'') 

for some constant L with < L < I. 

Proof: The first relation (|37|) follows from Lemma [3l Lemma |4] and Taylor's theorem. 

In order to prove (p8|) we first note that \x^\ < p — rj for i > 2 and < p + 77 (if 77 > is 
sufficiently small). Hence, by a second use of Taylor's theorem we get uniformly 

wu{x\u') = Ck{x'){u^ - l)y{x'f+^ + O (|u' - l\^y{\x'\f+^) 

and consequently 

i:k{x,u) = Y.\ Ckix^u^ - l)y{x^f+' + 0{\u- Ipydxp)^-) 

j>2 

= {u- l)Ck{x)y{x^f+^ + O {\u - lM\x\^f) . 
Here we have used the property that the sum 
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represents an analytic function in x and u, since due to |x| < p + e < 1 we have — > and 
therefore ^ impHes y{x^) ^ x*. FinaUy, since Ck{x) — C{x) + 0{L''), it also follows that Ckix) 
has a limit C{x) and the same order of convergence. □ 

With these auxiliary results we are able to get a precise result for Wk{x,u). 
Lemma 6. For u € Sfe and x € Ag (e, 9, fj as in Lemma\^ we have 

Ck{x)wo{x,u)y{x)'' 



Wkix, u) 



lCk{x)woix,u)'-^^ + 0{\u-l\) 



Proof: Since by Lemma[S]Sfc(a;, u) — 0{wk{x, u)), we observe that Wk{x, u) satisfies the recurrence 
relation (we omit the arguments now) 

Wk+i = y (e"'^+^'= - 1) 



y[wk + ^ + llk + 0{wi) + O(S^fe) 



= ynJk{l + ^ + (wl) + O (Efe)) (^1 + ^) 

Equivalently, we have 



Wfe+i WkWk+i Wk\ 2 

= — -]: + 0(wk) + o(— 
Wk 2 \Wk 

and consequently 



-1 2 V Wfe 



Wfe+l Wfe WkWk+l I \ Wk 

Thus we get by recurrence 

Wfc Wo ^ 21-j/ Vl-Ly V 1-2/^ 

Now we use again Lemma [5] to obtain 



Wo 



^ C,(x)Q+i(x)2/(x)^+i(l + 0(£|u- 1|)) 
= Cfc(a;)+0(?i-l) 

where Cfe(a;) denotes the sum in the penultimate line above. Observe, too, that wq ^^j^ 2 — Oil), 
if k\u — 1| < fj. Hence we obtain the representation 

(40) Wk = — k ■ 

l-Ck{x)-^^ + 0{\u~l\) 

Thus, it remains to verify that 1 — Ck{x) = 1/Ck{x). By using (l?f)) and ([5^ it follows that 

(41) Ckix) = Ec.(x^) = iC+ii^) ' Ck{x)) 

and consequently by telescoping 

fe-i 



. ^ Cj+ijx) - Cjjx) _ _L i_ _ _ 



C,{x)C,+^{x) Co{x) Ckix) Ckix) 
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Alternatively we can compare (^(11) with (I57| for m 1 which also shows 1 — Ck{x) — 1/Ck{x). 
This completes the proof of the lemma. □ 

The proof of Theorem |6] is now immediate. We substitute x = p (l + f^) (where s/n — >■ 0), 
u = e**/^ and set k — n^/n. We also use the local expansion y[x) = l — b^/py^l — x/ p+Odx — p\). 
That leads to 



y{x)'' = exp (-K,b^-ps) ( 1 + O 



\/n 

Finally, since the functions Ck{x) are continuous and uniformly convergent to C(x), they are also 
uniformly continuous and, thus, 

(42) Ck{x) ^ C{p) = h^p. 

Altogether this leads to 

wo{x,u)y{xY 



1 _ l-vi.)" ^0{\U- 1|) 

Ck(x) 2 l-v(x) ^1 1^ 

^ Ck{x){u^l){l~y{x))y{xf+^ 

1 - y{x) - Cfc(x)H«^(l - y{xY) + O (1^, _ i| . |i _ y^x)\) 
b'^p ity/^exp (— i/tfe^— ps) 



2%/"- -y^exp (1^6-^— ps) — sinh (1^6-^— pa ^ 



as proposed. 



5. The Finite Dimensional Limiting Distributions - Proof of Theorem H] 
For d = 1 (in Theorem |4]) we have 



Vn ^ 

(43) = / Vk (x, 

where the contour F = 7 U F' consists of a line 



1 /" ^ ^itl^^ '^^ 



7 = {x = p (^1 - ^^i^ j I -Clog^n < r < Clog^n} 

with an arbitrarily chosen fixed constants C > and cr > 0, and F' is a circular arc centered at 
the origin and making F a closed curve. 

The contribution of F' is exponentially small since for x G F' we have -^\x~^^'^\ — 

O ^n-^/^e" whereas j/fe ^x, e'*/^^ is bounded. 

If a; e 7, then the local expansion (IT4t is valid. Insertion into (|43l) . using ([5]), and taking the 
limit for n ^ co yields the characteristic function of the distribution of (i^'*) ^ desired. 

Now we proceed with d = 1. The computation of the two dimensional limiting distributions 
shows the general lines of the proof. An iterative use of the techniques will eventually prove 
Theorem [31 We confine ourselves with the presentation of the case d = 2. 

We have to show 
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Since the characteristic function of the two dimensional distribution satisfies 

^1 

Vn 



Vr, V / 



dx 

wn+l 



(45) 
where 

Wk,k+h{x,ux,U2) = ykM+h{x,ui,U2) - y{x) = Vk,h{x,ui,U2) - yix), 
we need to analyze the asymptotic behaviour of Wk^k+h{x, wi, U2) for k and h proportional to ^/n. 
Furthermore, note that y{x) and yk,h{x^ ui,ui) are analytic functions for a; € A and thus Wk.k+h is 

bounded for x € A. Hence the contribution of F' to the Cauchy integral (j45|) is O (^l/yn ■ e^'^'°s 

with some suitable constant c > and therefore negligibly small. Extending 7 to infinity, as it 
is required for the two-dimensional version of ()12p . again introduces an exponentially small and 
therefore negligible error. 

Thus it is sufficient to know the behaviour of Wk,k+h e'*^^ e'*^/ for a; G 7. 

Then, equation follows from the following proposition and ([S]). 

Proposition 1. Let K2 > ki > and ii,t2 &e given with |k2^i| < c anrf |K2i2| < c. Furthermore, 
set Ui — e'*'/ ^ and fc^ = [k^V^J « = 1, 2. Define s by x ^ p (\ + . Then we have 



(46) 



Wki,k2{x, Ul, U2) 



feV2^ 



2^/2 



2\/2 



2V2 



as n ^ 00, uniformly for x, ui and U2 such that ki\ui — l\ < c, k2\u2 — 1| < c, and 3 = (log^ n) 
such that 9fi(s) = —a where a > and sufficiently large but fixed. 

Proof: Note that Wfc^^fc^ (a;, 1, 1) = 0, 10^1,^2 (a;, "i, 1) = Wk^x^ui), and Wfe^^^^ (x, 1, U2) = 
Wk2{x,U2)- Thus, the first derivatives are given by 

-Wki,k2{x,ui,U2) =7fci(a;, 1), 



Wki,k2{x,Ui,U2) 

OU2 



U1—U2 — I 



lk2{x, 1). 



- U1—U2—I 

It is also possible to get bounds for the second derivatives of the form 0{k2y{x)''^), if x, -ui, and 
U2 are in the domain given in the assertion above. Hence, we can approximate Wki,k2{x, ui, U2) by 

(47) Wk,,k2{x.uuU2) - CkAx){ui - l)y{xp+^ + Ck2{x){u2 - l)y{xf'+^ 

+ 0{k2y{xp{\u^-l\'' + \u2-l\'')). 

Similarly (and even more easily, compare with the proof of Lemma [Sj we obtain a representation 
for 

Wkt,k2{x\u\,u\) 



Sfci,fe2(x,Ml,W2) = X! 



(48) 



= CkAx){ui - \)y{x^f'+' + Ck2{x){u2 - l)y{x^f^+' 
+ 0{y{\x\^p\u,^l\^+y{\xmu2-l\^)). 



In order to identify the asymptotic main term of Yiki^ix, ui, U2) note that k2~ ki ^ (^2 — Ki)y/n 
and |y(a;^)| < y(|xp) < 1 for |a;| < p + < 1. Moreover, observe that the terms ui — 1 and U2 — 1 
are proportional and hence 



(49) 



Ck2ix){u2 ~ l)y{xf- t2 I , , ^ 

CkAx){u,-i)y{x)^^ ^ (-('^^ - '^-)^^) 
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The Ki and the are fixed, so we can choose a such that the right-hand side in (1491) is different from 
1. This guarantees that the asymptotic main terms of (|Tf| . Cki{x){ui — l)y{x)''^ and Ck2{x){u2 — 
l)y(a;)'°^, do not cancel each other. 

By using the same reasoning as in the proof of Lemma |6] we get the representation 



WkiM 



1 - A, - ""■1"°^ \i; + Oi\u, -1\ + \u, 1| 



where 



fki = .fki{x,Ui,U2) 



fei-i 



Wo,k2-kiix,Ul,U2) ^ 



^iM-ki+iix, ui, U2) y{x) 



1=0 '^''-M-ki+iix, Ui,U2)wi+i^k2-ki+i+lix, Ui, -U2) 

Note that 

wo^k2-ki = uitk2-kt{x,U2) ~ y{x) 

= {ui - l)yix) + uiWk2~ki {x, U2) 
= U + W, 

where U and W abbreviate U — {ui — l)y{x) and W — wiWfc^-fcj (x, M2). The assumptions on Ui 
and ki given in the statement of Proposition [T] imply that by Lemma [5] we have (A >i B means 
that A and B have same order of magnitude) 

W X u^Ck2^k^{x){u2 - - o{u2 - 1) 

whereas C/xwi — lxu2 — 1. Thus we may safely assume that \W\ < ^\U\, so that there is no 
cancellation. 

Next by (gT]) we have 

Wi,k2-k,+e =Cdx){ui - l)y{xY+' + Ck2-k,+i{x){u2 - l)y(x)''-'''''+^ 

+ O ((fc2 - he)y{xY{\ui -1\^ + \U2- ip)) 

=Ci{x)y{xYU + Ck2-k,+e{x){u2 - l)y{xY--''''+' 

+ O {{k2 - kMxY{\ui -1\^ + \U2- Ip)) 

But since Ck2-ki+e{x) ^ Ck2-ki{x) by formula ([^5]) and Lemma [5] relates the second term to W , 
we obtain 

wi,k2~k^+i = {Ci{x)U + W) y{xY (l(fc2 - ki)\u2 - 1|)) . 
Hence, fk can be approximated by (for simplicity we omit the error terms) 

C,ix){y{x^)/y{x)Y+' 



fk,-UiU + W) ^ 
fci-i 

= {U + W)J2 



= iu + w) 
= 1 



{Ci{x)U + W){Ci+i{x)U + W) 

iCi+i{x)U + W)- {Ciix)U + W) 
{Cdx)U + W){Ce+i{x)U + W) 

1 1 



U + W CkU + W 

u + w 



CkU + W 

where we have used the formula (HT|) and telescoping. Consequently, it follows that 

{U + W)y''^ 
'^kiM ~ u+w i-v^ 

CkU+W 2 1-y 

r^m = iCkU + w)y^- 

2 l-y 
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,2 

Now we will approximate all the terms in ([50)) . First recall Ck{x) ~ ^ (formula (US)))- The 
assertion on a; = p (l + -^j implies a; — > p and hence — > 1. Thus we obtain 

The asymptotic expansion Q and Theorem [S] imply 

1 



y '1 - exp --Ki^V^ps 



W = MiWfe2-fci(a;, M2) 



2\/2 



it2by/p\ 
2V2 J' 



Applying all these approximations we finally get . 



□ 



6. Tightness - Proof of Theorem [5] 

In this section we will show that the sequence of random variables /„(t) = n~-^/^Ln(t^/n), t > 0, 
is tight in C[0, 00). By [3ll p. 63] it suffices to prove tightness for C[0, A]. Hence we consider L„(i) 
for < t < A^/n, where A > is an arbitrary real constant. 

By dJ Theorem 12.3] tightness of ln{t) = n-'^l'^Ln{ty/n), <t < A, follows from ti ghtness of 
J^ri(O) (which is trivial) and from the existence of a constant C > such that ([9]) holds for all 
non-negative integers n, r, h. 

The fourth moment in Equation (|9]) can be expressed as the coefficient of a suitable generating 
function. Indeed the generating function counting tree according to size as well as the quantity 
Ln{r) — Ln{r + h) is yr,h {x,u, ^) {yr,h{x,u,v) is defined by pTt ) since assigning a weight u to 
the vertices in level r and weight 1 /u to those in level r + h means that any tree having n vertices 
and with L„(r) — Ln{r + h) = to contributes to the coefficient of where m can be negative. 

The fourth moment can then be obtained by applying the operator [u-^j and setting u — 1 
afterwards. Therefore we have 



E(i„(r) -i„ 

Thus, (O is equivalent to 

d_ 

du 



Vn 



— +7 — 

du dv? 



VrJi 



X, u, — 
u 



u=l 



(51) 



7 



92 



6- 



g3 



9u4 



yr,h \X,U,- 



In order to prove (fSTj) we use a result from [18] saying that 



implies 



F{x)^0{{l-x/p) 
[x-]F{x) 



Oip 



where A is the region of (IT51) 

Hence, it is sufficient to show that 



(52) 



_9_ 

du 



7 



52 



6- 



93 



du^ du^ 



du^ 



yr,h X, U, ■ 



(xe A) 



O 



o 



for a; £ A and h>l. (Note that 6* < § implies that 1 
c > 0.) 

Now we define 



d^yk{x,u) 



and 



|y(a;)| > c-\/|l — x/ p\ for some constant 
d'yr,h {x,u, i) 



du3 



and derive the following upper bounds. 
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Lemma 7. We have 

/53) [1] <^ (1) uniformly for x € A, 

1 0{jx/p\^^ uniformly for \x\ <p 

and 

(54) = I O (^Tt) uniformly for x e A, 

[ 0(|a;/p|'") uniformly for \x\ < p, 
where L is constant with < L < 1. 

Proof: We already know that jI^\x) — Ck{x)y{x)^ , where Ck{x) — 0(1) and \y{x)\ < 1 for 
X €z A. Furthermore, by convexity we also have |y(a;)| < \x/p\ for < p. Hence, we obtain 
7i''W =0(|x/p|'=) for |a;|<p. 

The functions ^l^jiix) are given by the recurrence 

i>i 

with initial value 7g^jj(a;) — y{x) — 7/1 (x). To show this, just differentiate (|lip w.r.t. u and then 
plug in w = 1. Hence, the representation ^yf^ix) — ^\^\x) — ')'^^]^^{x) follows by induction. Since, 
-i\^\x) = (C(x) + 0{U))y{xY we thus get that 

7il(x)-ofsup|2/(xr(l-y(x)")|+L' 

However, it is an easy exercise to show that 

(55) sup|y(xr(l-y(x)'')|^of^y 

For this purpose observe that if a; G A then we either have \y{x) — 1| < 1 and |y(a;)| < 1, or 
< 1 — ?7 for some > 0. In the second case we surely have 

\y{xr {I - y{xf)\ < 2{l - ^0{U). 

For the first case we set y = 1 — pe"'' and observe that 

|l-(l-pe*'^)''| < + 

Hence, if r > 3/i we thus obtain that 

\y{xY{l - y{xf)\ < max (1 - pY ((1 + p)" - l) < < 

o<p<i ^ ' r — h r + h 

If r < 3/1 we obviously have 

|y(:.r(l-y(a:f)|<2<-i^ 



which completes the proof of ([55l) . Of course, we also have U = O (^tt^J- This completes the 
proof of the upper bound of "^^^\{x) for x e A. 

Finally, the upper bound l^^\{x) = O {\x/ pY) follows from ([5^ . □ 

Lemma 8. We have 

(56) ^[2] {x) = [ ^ [^^^ }) ^'^^f'^^^y f""" ^ S ^' 

I Oi^lx/pl^'j uniformly for \x\ < p ~ rj 



an 



d 



\2\ , X \ O ( min \ h, -. — } ] uniformly for x £ A, 

[ 0{\x/p\^) uniformly for \x\ < p — rj 

for every r] > 0. 
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Remark. By doing a more precise analysis similarly to Lemma [3] we can, for example, show that 
7^' (x) can be represented as 

fc 

(58) 7i''(:^^) = y(^)'E^M(^)y(2^r~\ 

where the functions Dk^i{x) are analytic in A. For every £ there is a limit Di{x) = limfc_j.oo Dk/{x) 
with 

DuAx)^D,{x)+0{L''+'), 
where < i < 1. Furthermore these limit functions Di{x) satisfy 

DiXx)^C{xf + 0{U). 
Since we will not make use of this precise representation we leave the details to the reader. 

Proof: The bound 7)?' (a;) — 0{\x/ p\^) (for \x\ < p — rj) and the bound "f^^\x) — 0{k) follow from 
Lemma |4l In order to complete the analysis for ^^^^ (x) we recall the recurrence derived from (jlOp 
(compare also with (|36|) ') 

(59) ^tUx)^y{x)Y,^^'^^\x-) + y{x) \Y.lt\xA + j;(a:) ^(z - 1)7™ (x^) 

i>l \i>l / i>2 

that we rewrite to 

(60) i^l,{x)^y{xW^\x)+hk{x), 
where 

(61) b,{x)^y{x)Y,z^^^\x^)+y{x) [y.1^^\x^)\ + y(x) ^(z - 1)7^(2:^). 

i>2 \i>l I i>2 

Note that for i > 2 we have |a;*| < p—rj. Therefore we can apply the second estimate of ([56]) (first 
and third sum of (|5T|)) and the estimate ([55)1 (the first for z = 1 in the second sum of (|5T|). the 
second for the other summands) and obtain then 

(62) bkix) = O (1) uniformly for x G A. 
Since 7o^'(x) — 0, the solution of the recurrence (|60)) can be written as 

(63) -f^^\x) = bk-iix) + y{x)hk-2{x) + ■■■+ y{xf-^bo{x). 
So (|62p implies finally 

'k-l 



7F'(a^) = oK:|y(x)P =0(^ 



\v{^)\ 



which completes the proof of (|56l) . 



The recurrence for (x) is similar to that of 7^,^' (x) 



(64) 7f|i, j^) = y{x) E + v^^) E + y(-^) E(* - 1)7S(^^) 



[2] I \ ^ ^ ■ [2] 

i>l \i>l I i>2 



\2] \2] 

with initial value 7q ^^(x) — 7^ (x). We again use induction. Assume that we already know that 
17^1(^)1 — Dr,h\x/p\'' for |x| < p — r] and for some constant I?r,/i- By (|57)) we can set Do,/i — Dh 
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which is bounded as /i — ?• oo. We also assume that hl'iix)] < C\x/p\'' for |x| < p~ rj. Then by 
we get 



Thus, we can set 



which shows that the constants Dr,h are uniformly bounded. Consequently 7|,^]j(a;) = O (|a;//9|'^) 
for |x| < p — rj. 

Next we start from ((64|) and assume that |7|.^]j(a;)| < £'r,h foi' a; e A. We already know that 
|7?I(a;)|<C^forxeA.Hence, 



i>2 



2 



Thus, we can set 

k + h 



with initial value I?o,h = = O [h) and obtain a uniform upper bound of the form 

Dr,h ^0{h). 

[21 

Consequently 7^ ^(x) = O [h) for a; S A. 

Thus, in order to complete the proof of (|57|) it remains to prove l^^\{x) = O (1/(1 — |y(a;)|)) for 
X e A. Analogously to the way we obtained from ([5^ we obtain from the representation 

of l^^l{x) as 

(65) lf}hi^) = 7o^i,(2^) + Ck-i,h{x) + y(a;)cfc_2,h(a;) H h y(a;)''"^co,/i(a:), 

where 

c,,.(x)=y(x)^»7g(x') + y(x) (^^^(^Ol +y(^)E(*-l)^S(^')- 

i>2 \i>l / i>2 

Observe that there exists 77 > such that |a;*| < p — rj for i > 2 and a; G A. Hence it follows 
in a similar fashion as we showed (15^ that Cj,;i(a;) = O (1) for a; e A. Since 7^ \_^{x) = 7)^ '(x) = 
0(1/(1 — |j/(a;)|)), we consequently get 

1 



7|l(^)=7i''(-)+o(^ 



\y{x)\J Vi-|y(^)l 



□ 
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Remark. Note that the theorem our proof rehes on, namely [5J Theorem 12.3], actually requires 
the existence of a > and /? > 1 such that 

E |i„(r) - L„(r + ^01" - O {{hV^f) . 

The estimates of Lemma [5] already prove that 

E (L„(r) - L„(r + h)f - O {h^ . 

Unfortunately this estimate is slightly too weak to prove tightness. Third moments are technically 
unpleasant they attain positive and negative signs. So we actually have to deal with 4-th moments. 

Before we start with bounds for ^'^^\x) and l^^\x) we need an auxiliary bound. 

Lemma 9. We have uniformly for a; e A 

(66) El^?i(^hS(^)|-0(/.^). 



Proof: We use the representation (|65p . where we can approximate Cj^h{x) by 



(r + hf 



uniformly for x £ IS. with some constant L that satisfies < L < 1. Furthermore, we use the 
approximation 

7S(a;) = C{x)y{xni - y{xf) + O [U) 
that is uniform for x G A. For example, this shows 



r>0 



|y(^)l 



Now observe that for a; e A there exists a constant c > with |1 — y{x)\ > c(l — |y(a;)|). Hence it 
follows that 



|1- 2/(3^)1 



O 



1 - y(x)'' 



1 - y{x) 



0{h) 



and consequently 



Similarly we get 



r>0 



r>l 



Y,y{xy-^~'c,,u^ 



j<r 



< 



j>0 



r>j 



E \cjAx)\ \y{x)r-' (\C{x)\ \yix)\'^+^\Ly^ + o {\y{xW+'L^) 



yj>0 

= O (h') . 



ij + hf 



Hence, we finally obtain 

E l7S(^)7S(2:)| < E iT^l^;)! + E 17^(^)1 



17;. 



Y,y{xr~=-'cM 

<r 

Y^y{xr-^-\,,u{x) 



o{e) 
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where we used the fact |7jj {x)\ — O (h) from Lemma HI □ 
Lemma 10. We have 



(67) il: 

and 



[3]^^^ ^ I O (^min{fc2, __^|^}^ uniformly for x € 



O uniformly for \x\ < p — r] 

O (min{h'^, i) uniformly for x € A, 

O (}x/ p\^) uniformly for \x\ < p — rj 



for every 77 > 0. 

Proof: The recurrence for 7]^ (x) is again obtain by differentiation of (jlOp and given by 

i>i yi>j y yj>i y \i>i 

(69) +3y(x) (E^Fk^^) I (E(*-l)7i^'(^^) I +3y(a;)E*(*-l)7i''(^^) 

\i>l J \i>l 

+ y{x)J2i^-m-2h]^\, 



, j>l / \i>l / i>l 

ix') 



By inspecting the proof of Lemmas [7] and [8] one expects that the only important part of this 
recurrence if given by 

(70) (x) = 2/(x)7f 1 (x) + yixWk^ {xf + 3y{xW^^ {x^^^ (x) + R, 

and Rk coUccts the less important remainder terms that only contributes exponentially small 
terms. Thus, in order to shorten our presentation we will only focus on these terms. In particular 
it is easy to show the bound "f^^\x) = O {\x/p\'') for \x\ < p — rj. (We omit the details.) 

Next, since y{x)j^^\x)'^ + 3y{x)j^^\x)^^^\x) + Rk — O (k) for x G A, it directly follows that 
7Fl(:^)=0(fc2). 

Now we proceed by induction and observe that a bound of the form |7[^'(a;)| < Ek/{\ — \y{x)\) 
leads to 



and consequently to E^+i < Ek + O (1). Hence, E^ ^ O (k) and J^^^) = ^ ('^/(l ~ Ivi^)])) 
Similarly, the leading part of the recurrence for 7^ ^(x) is given by 

(71) i^rlij^) = yi^hZi^) + yi^hlii^f + M^Wrli^h^li^) + Kn 

where 

drA^) = 2/(x)7i^I(x)3 + iy{x)^^^l{x)^'^^l{x) + Rr^n - O (h) 
and the initial value is given by 

7a(x) = -7i''(x) - 37r'(2:) = O (^min !^h^ ^3]^}) ' 
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Note that we also assume that J^^]^{x) — O {\x/p\^) for < p — r] (which can be easily proved). 
Consequently it directly follows that 



^^rhi^) = "f^oh(^) + dr-i.h{x) + y{x)dr-i,h{x) H h y{xY ^do,h{x) 

= 0^ ^ 



.^~\y{x)\. 

Next observe that Lemmas [7H9] ensure that 

Y,\d,,h{^)\^0{h') 

uniformly for a: € A. Hence, we finally get 

which completes the proof of Lemma (TU] □ 
Lemma 11. We have 

(72, = { «(t^) a, 

[ O (|a;/p|'^) uniformly for \x\ < p — r] 

and 

,73, .;i«J''(t^) /-"A. 

I 0{\x/p\^) uniformly for \x\ < p — rj 

for every rj > 0. 

Proof: The proof is very similar to that of Lemma [TOl First, the recurrence for 7^,^' (x) is essentially 
of the form 

(74) l^lA^) = yixW^\x)+yixW^Hx)^+4yix)^'i\xW'\x) 

+ 6y{xW^\xfj^^\x) + 3y{xh^\xf + i?,, 

where Rk collects all exponentially small summands. We assume that we have already proved the 
upper bound "f^^H^) = O (|a;/p|'^) for |x| < p—rj. Now, by induction and the assumption |7[!^'(x)| < 
Fk/{1 — \y{x)\) and the known estimates jI^\x) = O (1), ^]^Hx) = O (min{fc, 1/(1 — |y(a;)|)}), and 
7['^(x)=0(fc/(l-|2/(x)|)) weget 

hll. wi + + (^) + (^) + m 

and consequently Fk = O {k'^)- 

Finally, the essential part of the recurrence for ^y^^j^ix) is given by 

(75) 7^1, J^) = y(x)7?i(x) + Vi^hllixr + 42/(x)7^i(:^)7S(x) 

+ 62/(.T)7?I(x)^S(x) + 3y{xWrlixr + i?.,, 

where 

e,,,(x) = y(:r)7W (:.)^ +4y(x)7W (x)7S(2;) 

+ 6y{xWrlixf7^'lix) + SyixW^lix)' + Rr,H. 
As above, Rr^h collects all exponentially small terms. Thus, 

lt\i^) = 7o'^i(a;) + er-i,h{x) + y{x)er-i,h{x) H h y{xY^^ eo^h{x). 



THE SHAPE OF UNLABELED ROOTED RANDOM TREES 



25 



If we use the known estimates "fl''\{x) — O (1), J^rhi-"^) ~ ^ (^)' ^^^^ "f^rhi-"^) ~ ^ (^^) which gives 
dr,h = O (ft-^) and the initial condition 



we obtain 



This completes the proof of Lemma [TT] □ 
The proof of ((52|) is now immediate. As already noted this implies ([5T|) and proves Theorem [51 

7. The Height - Proof of Theorem [2] 

(k) 

Let j/n denote the number of trees with n nodes and height less than k. Then the generating 
function ykix) — J2n>i Vn^x" satisfies the recurrence relation 

yo{x) = 

yk+i{x) = xexp I ^ ^(^iifj. j ^ fc > 0. 

Obviously yk{x) = yk{x,0) where the function on the right-hand side is the generating function 
of (jlOp which we used to analyze the profile in the previous sections. So Wk and Sfc could be 
defined accordingly. However, the proof given here relies heavily on the seminal work of Flajolet 
and Odlyzko [20] on the height of binary trees. Therefore, to be in accordance with the notation 
used there, we work with the opposite sign and set 

ek{x) = y{x) - yk{x), 

that is ek{x) — —Wk{x,0). Then ek satisfies the recurrence 

(76) ek+i{x) = y{x) (l - e-^'^^^)-^'"'-)) , 

where 

efe(a;*) 

i>2 



(77) Ek{x)=y2^^^~Ek{x,0). 



The function ek{x) is the generating function for the number of trees with height at least k. 

The proof of Theorem [5] follows the same principles as the proof of the corresponding properties 
of the height of Galton- Watson trees (see [101 [IS]) - However, the term Ek{x) needs some additional 
considerations. 

Proposition 2. If ek{x) satisfies 
(78) ek{x) = 



y{xY 



1 l-y{x)>' 

2 l-y{x) 

for X € Ag, then Theorem\^ follows. 



o(min{logfc,log^3^} 



Proof: Since ek{x) is bounded in A \ Ag, only the local behaviour near x — p determines the 
asymptotic height of Polya trees. The shape (|75)) of Ckix) precisely matches that of the corre- 
sponding quantity for simply generated trees. Flajolet and Odlyzko showed that ([78]) implies ([6]), 
see [101 P- 204] where this argument was used to derive the average height as well as the other 
moments of the height of simply generated trees. □ 



26 



MICHAEL DRMOTA AND BERNHARD GITTENBERGER 



Proposition 3. Suppose that for x e the estimate 

(79) Wh-Oi^L^) 
holds for some L < 1 and that 

(80) \ek{x)y{x)''\^0{l/k). 
Then we have (j78p for x € . 



Proof: Equation (j76p can be rewritten to (omitting tire argument x) 

ek+i = y^k\i~ Y + [^k + — 



rasp, to 



This leads to the representation 

Recall that eo = y{x). By ([79]) and = O (1) (for x £ Ag), a consequence of LemmalU this 

implies 

(82) efe 



2 1-y 

and dSni) yields ([751) after all. □ 

Remark. Note that the proof above does not make explicit use of the domain of x. Thus the 
implication of Proposition [3] is still true if we write, for instance, < a; < p instead of x G A^ in 
([75)1 . ([7^ . and ([5D|) . We remark that we will use such modifications of Proposition [3] in the sequel, 
though we do not state several almost identical propositions differing only in the domain of x. 

Note that (ISTt and ([75]) can be made more precise. Set 

el {e-^" - 1) 



(e^^fc - 1) (1 - e-<^>'-E'^y 

and define a function h(v) by 

1 + 7: + v'^h[v). 



V 



1 - e^" 2 
Then the recurrence Ck+i = 2/(1 — e^'^k-Ek^ rewrites to 

y 1 1 1 I \ 

= ^ o + ekh[ek) + 

Cfc+i e/c 2 e^. 

and leads to the explicit representations 
(83) 

£<k i<k ^ 

and 

(84) efc - ^'^ 



i + i^ + E.<.e.Me.)2/^ + E.<.§?/^ 



This formula is a refinement of ([78]) since it makes the error term explicit. We will use it in the 
sequel. 

Furthermore, note that if we just assume Cfe — )■ and = o(efe) as A: — )■ 00, then 
(85) Sk ~ -Ek and h{ek) = 0{\). 
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We start our precise analysis with an a priori bound for efe(x). The next step is proving (j78p 
for < X < p. Then we will, little by little, enlarge the allowed domain for x and arrive finally at 
Proposition [3] as stated above. 

Lemma 12. Let \x\ < p. Then there is a C > such that 



\ek{x)\ < -j= 



C - 



Proof: Obviously, we have 

|efe(x)| = ^(y„-yW)N"<^y„ 



. - x\ 



The assertion follows now from ?/„ ~ cp~^n~^/'^ for some constant c > 0. □ 

Lemma [T^l applies to Ek{x). 
Corollary 3. Suppose that \x\ < Then there exists a constant Co > with 



\Ekix)\ < ^ 



„2 



fc 



X 
P 

Remark. Observe that in the definition of Ek{x), Equation (|77p . the arguments in the sum are 
raised to a power of at least 2. Therefore Ek{x) is an analytic function in the domain \x\ < ^/p 
and not only in the smaller domain < p. 

The next lemma shows that ek{x) behaves as expected if x is on the positive real axis. 

Lemma 13. Suppose that < x < p is real. Then |7i$|j holds in this domain for x. 



Remark. As remarked before, we will show a weaker version of (|78p. During the proof we will make 
use of the analogous weaker versions of ([79]) and ((80|) . Therefore, in the following proof a reference 
to one of these formulas means the referred formula, but for < a; < p and not for x € A^. 

Proof: Let ek{x) be defined by eo(a;) — y{x) and by Ck+iix) ~ y{x){l — e"'^'''^^^) (for k > 0). Then 
ek(x) is precisely the analogue of for Cayley trees, a class of simply generated trees (preceisely: 
the class of labelled rooted trees). So ek{x) behaves like (f78|) in A. 

However, if < a; < p then we obtain by induction that ek{x) > ek{x). Hence, by combining 
(|78p with the upper bound from Lemma [T^] we have 

Ek{x) ^ Ek{x) ^ ^ 
ekixy ~ Ckix)'^ 

for some L with < L < 1. Thus is satisfied. 

In order to show the second assumption of Proposition |3] note that by Lemma [T^ ek(x) 
is even exponentially small for x < p. For the case x — p observe that (|79p in conjunction with 
Lemma fT2l guarantee (|85|) . Applying this to (|84t implies 

y{xf 



ek{x) 



This equation yields 

(86) ek{p) ^ I 

and completes the proof. □ 

The analysis of ek{x) for complex x with |a;| < p is not too difficult. The next two lemmas 
consider the case |a;| < p and |x — p| < e and the case |a;| < p — e. 

Lemma 14. There exists e > such that \78^ holds for all x with \x\ < p and \x — p\ < e. 
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Proof: First recall that |efe(a;)| < C/Vk and \y{x)\ < 1 for x 7^ p. Moreover, in the proof of the 
previous lemma we showed ek{p) ^ 2/k. Hence ([50)) is satisfied. 

Suppose that we can show that \Ek/e1\ < 1 or, equivalently, \Ek/ek\ < \ek\ < C/Vk. Then it 
follows that 

1 - e-'"'-'^'' 



|efe+i| = \y\ \ek\ 



= \y\\ek\ (1 + ^ + 
> |y||efc|l-Ci|efe| 



{ek + EkY' 
efc 



(87) 



> \y\ |efc|e 



where Ci,C2 are suitable constants. 

Now we choose fco sufficiently large such that 



k 

hold for all k > kc,. By continuity, (|86|) implies the existence of an e > with |efeo(a;)| > ^ and 
> p^^^ for < p and |a; — p| < e. These assumptions imply 

\eko{x)\ > ^ > e-2C^^ > |y(a;)|'^-°e-2C2V^ 
and by Corollary |3] (since |a;| < p < this is applicable) 

The goal is to show by induction that for k > ko and for |x| < P and |x — p| < e 

Ek 



< 1. 



Assume that 
implies 



lefel > Ij/l'^e-^c^^^ and 



is satisfied for k = kg. Now suppose that (|88p holds for some k > ko- Then ([87 

\ek+i\>\y\\ek\e-''^'"'' 
> |y|'=+ie-3C2/fe+T. 



Furthermore 



\Ek. 



+11 



<CoP^-+'|?/| 



fc+l L,|-2fc-2„4C2VFFT 



Hence, we have proved ([88[) for all fc > /co- 
in the last step of the induction proof we also obtained the upper bound 



\Ek\ 



< CoP^^^e^^-^ 



which is sufficient to obtain the asymptotic representation ([78p . 

Lemma 15. Suppose that \x\ < p — e for some e > 0. Then we have uniformly 
(89) ek{x) = Ck{x)y{x)'' = {C{x) + o{l))y{x)'' 



□ 



for some analytic function C{x). Consequently we have uniformly for \x\ < p — e 
(90) Ek{x) = Ck{x)y{xY = {C{x) + o{l))y{x^)'' 

with an analytic function C{x). 
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Proof: li \x\ < p — e then by Lemma [T^ we have |efe(a;)| < ek{p — e) = O ^^1 ^ j;^ ^ • Thus, we 

can replace the upper bound |efc(x)| < C/Vk in the proof of Lemma [14] by an exponential bound 
which leads to a lower bound for efc(x) of the form 

\ek{x)\ > co|y(x)|^ 

Hence, by using (IM)) the result follows with straightforward calculations. 

In order to show the second assertion we start from (|77p and insert (|89p . Then we obtain 

(91) /£|^^^ (a(..) + .(i)M..)'- ^ 

y i>3 ' J 

Since C{x) is analytic, it is bounded in the compact interval [0, p — e]. Furthermore, observe that 
using the bound from Lemma [1] it is easy to see that the series in (|9ip is uniformly convergent. 
Hence the representation ([TO]) follows. □ 

The disadvantage of the previous two lemmas is that they only work for |a;| < p. In order to 
obtain some progress for |a;| > p fix a constant C > such that 



~ Vk\ P 



for all fc > 1 and for all |a;| < y/p. 

Lemma 16. Let a; G A and suppose that there exist real numbers Di and D2 with < Z^i, D2 < 1 
and some integer K > 1 with 

(92) \eKix)\<D^, \yix)\^^^-^ < D2, D^D2 + e^^-^(^-^) < D,. 

Di yjK \ P J 

Then we have \ek{x)\ < Di for all k > K and 

efc(x) = 0{y{xt) 
as — >■ cx), where the implicit constant might depend on x. 

Remark. Note that the second inequality in (j92p implies |y(a;)| < 1 and thus an exponential decay 
of ek{x). Hence, in particular, the assumptions of this lemma imply that (jSOp holds. 

Proof: By definition we have ck+i — y{l — e^'^'^^^^). Hence, if we write e^^'' ~ l + Rk we obtain 

\eK+i\ < IcKWyl^-^^ + e^'-^RK. 
\eK\ 



If ((92|) is satisfied then it follows that 



VK \ p J 

Now we can proceed by induction and obtain \ek\ < Di for all k > K. Note that D2 < 1 and since 
x S A we have |a;| < p + e, Corollary [3] implies 



C /IxP^"^ 



\Ek{x)\ < ^ 



P 



Moreover, Rk = e^^" - 1 = O (|£;fc|) and hence J2k = 0{1). Thus we have 

Ck = 0{D'^). 

If we set ak = y{l — e^'^'')/ek and bk ~ —ye^'^''Rk we obtain the recurrence 

efc+i = efcflfc + bk 
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with an explicit solution of the form 

K<i<k K<j<k j<i<k 

Since — ©(Dj), we have 

n a, = 0(2;(xf-^) 

j<i<k 

and hence by Corollary |3] 

n = O = O {y{xfL^) 

j<i<k 

for some L with < L < 1. Hence, 

(93) ek^O {cKvix)''-'') + O {y{xfL^) = O {y{xf) 



□ 



Recall that ek{x) — > for |a;| < p. Using Lemma Il6l we deduce that — > in a certain region 
that extends the circle |x| < p. 

Lemma 17. For every e > there exists ^ > such that ek{x) at an exponential rate, if 
\x\ < p + 5 and |a; — p| > s. 

Proof: If we show that Lemma [16] is applicable, then we are done. First observe that since the 
assumption on x implies < 1, for sufficiently small Di there exist Z?2 < 1 and K (sufficiently 

large) such that the second and third inequality if (|92|) hold. Next note that |ei<-(x)| < exip) for 
|a;| < p and thus |e/f(a;)| < 2eK{p) by continuity. By (|86|) we have exip) ^ ^ we can make |ex(a;)| 
arbitrarily small and the proof is complete. □ 

Now we turn to the most important range, namely for x e A and |a; ^ p\ < e. 

Lemma 18. There exists £ > and a constant ci > such that for all x G the conditions 
\92\) are satisfied for 

Uarg(y(x))L 

and properly chosen real numbers Di,D2. Consequently, ek{x) at exponential rate. 

Proof: Suppose that arg(x) and a.Tg{y{x)) are positive and that K{x) — 
[ci/l arg(y(a:;))|J , where ci = arccos(l/4) — £i and ei > is arbitrarily small. Note that 
K{x) can be made as large as we desire, since y{x) ^ 1 and therefore arg(y(a;)) is small for small 
e. 

Now fix an integer kg and e small enough to guarantee fcp < K(x). Moreover, fix two small 
positive real numbers £2 and £3. First we will prove by induction that for kQ < k < K{x) we have 
(for sufficiently small £) 

(94) |efc(a:)|<£2, 

(95) arg(efe(a;)) < k arg(?/(x)) + £3, 

(96) |efc(a;)| < c '^\^' for some c> 0. 

k 

The first step of the induction proof is to show (|M1) . Observe that due to the choice of ci and 
formula (|95|) of the induction hypothesis we have 

(97) < arg(efe(x)) < arccos(l/4). 

This implies |efc+i| < |efc| + |i?fe| and consequently, by the second statement of Lemma [TSl and the 
property ek„{p) ^ c/ko (compare with (j86p ) it follows that 

(98) \ek{x)\ < \eko\ + l^il < £2 

ko<i<k 
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provided that fcp is chosen sufBciently large. 

Next we show (PSI) . We start with ([75]) and obtain 



(99) 



1 



efc 



efe 



Note that by (I94p the first of the two error terms is much smaller than Cfe. In order to estimate 
the second error term, note that by the second statement of Lemma [TSlwc know Ek{x) = {C{x) + 
o{l))y{x^)^ . Combining this with (|94| we obtain 



O 



ky{x 



2\k 



0{L^) 



with some < i < 1. Since \y{x'^)/y{x)\ < p (& consequence of the convexity of y{x) on the 
positive real line) whereas \y{x)\ ~ p, we can have L < \y{x) \ provided that e is small enough. But 
this together with ([M)) implies that also the second error term in (IM)) is small in comparison to 
Cfc. Hence (|97)) implies that the argument of the last factor in (l99l) is negative. Thus we conclude 



arg(efe+i(a;)) < arg(y(a;)) +arg(efe) 

< (fc - /co + 1) arg(?/(a:)) + arg(efeo(a;)) 

< fcarg(y(a;)) +£3 

c/ko (compare with 



and continuity 



where the last inequality follows the fact that by Cfe^ (p) 
we can always achieve | arg(efeo(a;))| < £3. 

The third step is to prove the lower bound (|96l) for ek{x) for fco < k < K{x). of the form 
|efe(a;)| > c\y{x)\'^ /k for some c > 0. By hemn\a,\Y^Ehix) — {C{x) + o{l))y{x^)'^ behaves nicely, if 
\x — p\ < £. Suppose that \x\ > p and x e A^. Since arg(y(a;^)) is of order arg(?;(a;))^ we deduce that 
arg{Ek{x)) = 0{aig{y{x))) for ko < k < K{x). In particular, it follows that (for ko < k < K{x)) 

\ek+i{x)\ - \yix)il - e--a-)-^^'=(-))| > |y(2;)|(l - e-l'^^(-)l). 

Treating the nonlinear recurrence ak+i = 1 — e^°'= with the methods of de Bruijn [7l p. 156], it is 
possible to show inductively that afe ~ f and thus |efc(x)| > c\y{x)\''/k for some c > 0. 

The last task is to find Di and D2 such that the conditions ([M)) are satisfied for k — K{x). In 
order to do this, we first show that in the formula 



(100) 



..k—kc) 



efe 



1 



1 i-y 
2 



Eko<i<k ^Mei)y^ ■ 



where fco is fixed, the second term in the denominator dominates. Since fco is fixed, the first term 
is bounded. For estimating the third term note that by (|M)) the terms in the sum can be made 
arbitrarily small if fco is chosen sufficiently large. Finally, due to the already obtained bounds 
|efe(a;)| > c\y{x)\'^ /k, Ek{x) — 0{y{x^Y), and the property Si ^ —Eg the last term satisfies 

Si 
e? 



ko<l<k 

and therefore it does not contribute to the main term, either. Summing up we have 

|efe|< J.^' , (1+^4) 



2 l-y 



for an arbitrarily small £4 > 0. 

We set arg(p — x) — 9, where we assume that 9 G [— f — £5, -I + £5] (for some £5 > that has 
to be sufficiently small), and r = b\p — a;|^/^, where b is the constant appearing in (g]). Then we 
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have 

Q 

\y\ = \ -rco&- + 0{r^), 

Q 

log|j/| = -rcos- + 0(r2), 

Q 

arg(?;) = -rsin- + 0{r'^). 
Hence with k = K{x) ~ [ci/| arg(2/)|J we have 

|y'=-'=0 I ^ g-ci cot(e/2)+0(r2) < g-ci cot(f + ^)+0(r^) < g-ci (^-^ _ ^^-j 

for some arbitrarily small Eq > (depending on £5). Consequently 

g— arccos(l / 4) 
l^fcl < Dl ■■= 2 ^_^_arccos(l/4) Kl + ^t) - cV, 

where ey > can be chosen arbitrarily small. Moreover 

|y| = 1 - r cos ^ + 0(r2) < 1 - -^^(i _ g^) 

for some (small) eg > and consequently 

g-Di _ 1 / 2 e~ arccos(l/4) \ 

Thus, we are led to set 

1 / 1 g- arccos(l/4) \ 

^2 1 - t: ^ 7T7TT \ r = l- d'r 

and the first two conditions of are satisfied if r is sufficiently small. Since Di ~ D1D2 = c'c"r^ 
we just have to check whether 

C 



X 

P 

However, since k = K[x) > Ci^/2r~^ the left hand side of this inequality is definitely smaller than 
c'c"r^ if r is sufficiently small. Hence all conditions of ([M)) are satisfied for k = K{x). □ 

Lemma 19. There exists e > such that ^78^ holds for all x with x € with \x\ > p. 

Proof: We recall that the properties Ck = 0{l/k) and ([75]) imply ([75]) . By Lemma [T51 we already 
know that the first condition holds. Furthermore, we have upper bounds for Ek (see Lemma . 
Hence, it remains to provide proper lower bounds for Cfe. 

Since we already know that Cfc — J' and \Ek\ — 0{L'') (for some L < 1), the recurrence ([7^1) 
implies 

\ek+i\>{l~S){\ek\^\Ek\) 

for some 6 > provided that a; £ A and — p| < e. Without loss of generality we can assume 
that L < {l-6f. Hence 

lefcl > (1 - 5)'= - 5] \E,\{1 - 5f-' > co(l - 5f 
e<k 

for some constant cq > 0. Consequently 

o 



Ek 



As noted above, this upper bound is sufficient to deduce ([75| . □ 
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